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ABSTRACT 


An analysis of the fully coupled viscous, radiating flow past an ablating 
blunt body at hyperbolic entry conditions is presented. A detailed thermodynamics 
computation, as well as a realistic radiation transport model, is included in 
this analysis. A locally nonsimilar approach is employed to solve the conservation 
equations away from the stagnation point. The validity of the locally nonsimilar 
approach is demonstrated for some nonablating cases. Sample calculations are 
made for the typical flight condition of a Jovian entry probe. The effects of 
the downstream injection of the ablation products of a carbon heat shield 
on the flux distribution around the body are discussed in detail. It is found 
that most of the radiative energy is absorbed by the injected carbon gas and. 
dumped into the wake. A companion report, (NASA CR-2235) is available of the 
SL-4 Code and instructions for operation. 



NOMENCLATURE 


B(y, v) Planck . function 

C Mass fraction of the ablation gas 


D s (°) 


P C . 

1 

Effective binary diffusion coefficient 

Binary diffusion coefficient behind the shock at stagnation point 

u/u , normalized tangential velocity 
s 




Radiative flux 


Radiative flux at wall 


F (M=0) 
rw 


Radiative flux at wall when injection rate is zero 

Total enthalpy, also normalized total enthalpy, normalized by H 


h Static enthalpy, also normalized static enthalpy, normalized by H 

s 

l(\Uy,Q) Spectral radiation intensity 
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Shock curvature 
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00 
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A 
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H 


Total injected mass flow per unit area as a function of r 

Pressure 

Prandtl number 

Radiative flux divergence 

Distance measured from axis of symmetry 

Shock radius of curvature at stagnation point 

Reynolds number P co u ao R s /d s (°) 

Schmidt number ^(oJ/p^D (o) 

Temperature 

Temperature behind the shock at stagnation point 
Tangential velocity, velocity component parallel to the shock 

Free stream velocity 

Velocity component perpendicular to the shock 

Distance parallel to the shock measured from stagnation line 

Distance normal to the shock 

Distance along the line of symmetry from the shock stagnation point 
Viscosity 

Viscosity behind the shock at stagnation point 
Absorption coefficient 
Elemental gas mass fraction 
Parameters defined by Eqs. (ll) and (12) 

Shock layer thickness 

Transformed variable § = x 

Normalized stream function defined by Eq. ( 7 ) 


v 



♦ 

Stream function 


€ 

Stagnation point density ratio across the shock 

V 

Frequency 


P 

Density 


CO 

Parametric variable defined by Eq. 

(21) 

Subscripts 

s 

Quantities at shock 


b 

Quantities at body 


w 

Quantities at wall 


00 

Freestream conditions 


? 

Partial derivative with respect to 

variable g 


vi 



1.0 


INTRODUCTION 


The purpose of the present study is to investigate the effect of surface 
blowing on the thermal environment of a vehicle entering planetary atmospheres 
at hyperbolic speed. Most of the earlier studies which included the effect 
of ablation products are confined to the stagnation region (Refs. 1-5)- The 
only solutions obtained away from the stagnation region appear to be that of 
Olsted's (Ref. 6). Olsted adapted a two layer inviscid model and treated air 
to air surface injection only. Since no diffusion mechanism is allowed in 
the inviscid model, it probably is too crude an approximation for treating the 
effect of ablation products. In the present report, we study the thin-radiating 
shock layer about blunt bodies with surface mass injection, and do not restrict 
the analysis to the stagnation region. Viscous effects, as well as diffusion 
of the injected materials, will be considered. 

In our previous work, reported in Ref. 'J, solutions for a thin radiating 
shock layer about blunt bodies without blowing were obtained by a locally 
nonsimilar approach. Perfect gas and constant transport properties were 
assumed and only the continuum radiation was considered. Due to these 
simplifications, especially the nonblowing restriction, the results were not 
quantitatively realistic. Validity of the locally nonsimilar approach in 
solving the thin shock layer problem was sufficiently demonstrated, however. 

The work reported herein is a continuation of that of Ref. r J. The additional 
tasks were to incorporate a realistic radiation transport model as well as more 
detailed thermodynamic properties; to allow the ablation products to diffuse 
into the shock layer and to rederive the governing equations so that mass 
injection rates variable with distance around the body could be treated. 

The approximate but still complete treatment of the radiative transport 
in mixtures of hydrogen, helium, carbon, nitrogen and oxygen developed by 
Wilson (Ref. 5) was adapted to the present work. In this treatment, continuum 
radiation is approximated by a twenty-six-band model and atomic line radiation 
is represented by eleven- line groups. The radiation field, i.e., the flux 
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and its divergence, is determined by the solution to the one -dimensional 
radiative transfer equation. The validity of using one -dimensional radiative 
transfer equations for thin shock layer problems has been demonstrated as 
adequate in Ref. 8. 

The shock layer is assumed to be in thermochemical equilibrium. 

Hence, given pressure, temperature (or enthalpy) and elemental species concen- 
tration at any point in the shock layer, other thermodynamic variables can be 
completely determined by a Free Energy Minimization procedure. This is imple- 
mented in the present work by the FEMP Computer Program, described in detail 
in Ref. 9- Transport properties for a gas mixture, such as viscosity and heat 
conductivity, are computed by the approximate formulas of Refs. 10 and 11. 

The method of solution employed here is an inverse one, that is to say, 
the shock is completely specified and the body shape is a part of the solution. 
For a specific body shape, therefore, iterations would have to be performed on 
the given shock shape. Fortunately, the problems we will be interested in are 
thin shock layer problems; namely, the shock shape does not deviate much from 
the body shape. A good first guess on shock shape can be made and iteration 
largely can be avoided. 

Calculations for a body shape and flight condition, which represents a 
typical entry of a Jupiter probe, were made. The injected gas is assumed to be 
pure carbon, and the injection rate is a monotonically decreasing function of 
shock radius. The radiative heating distribution around the body is presented. 
This distribution is compared with that of no surface mass injection. 
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2.0 


MATHEMATICAL DEVELOPMENT 


2.1 Governing Equations 


The conservation equations for an axisymmetric thin shock layer, 
written in a shock oriented coordinate system (Fig. l) are: 


Continuity 


Spur 

dx 


= 0 

Sy 


(i) 


x- momentum 


d / du N du ~ du ^ . 

P ^ d^ } " PU d^ ' pVk d^ + pkUV 


i *° (2) 


y-momentum 


specie ^ 


sr)- 


. 1 X du‘ dH 

L_ Pr dy ^ " pU dx 


f 

d r u dH ti , 1 \ du" 

energy 5 y * 2 ( 57 

aC v - ”v 

kpT 5 T ' ° 


dH r dH 

55 * kpv + 


dC 


PU dT 


(3) 

( 4 ) 

(5) 


The y-momentum equation, Eq. ( 3 ), has been simplified according to Maslen's 
approximation (Ref. 12). The validity of Maslen's approximation has been 
examined in Ref. 7- In the energy equation, Eq. (4), the Lewis number is 
assumed to be unity. The symbol Q represents the radiative flux divergence, 
the precise form of Q will be discussed in the next section. 

In the specie conservation equation, Eq. ( 5 ) , C v denotes the mass fraction 
of the ablation gas. Here we have assumed a binary diffusion process, i.e., 
the ablation gas diffuses as a whole into the shock heated atmospheric gas. 

The detailed derivation of Eq. ( 5 ) can be found in Ref. 1 or Ref. 4. This 
binary diffusion model will be valid if the gas mixture can be divided into two 
groups of species, each with about the same atomic or molecular weight and 
about the same mutual collision cross section. Otherwise, multicomponent 
diffusion should be considered. For the sake of simplicity, the binary model 
was used in this study. 
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Fig. 1. Sketch of Geometry 
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Once the mass fraction for ablation gas is known, the mass fraction 
for elemental gas i is given by 


a. 

1 


(o' ).C + (or ).(1-C ) 

vi v a i v 


( 6 ) 


Where (a^) . is the mass fraction of element i in the ablation gas, and 

(or ) is the mass fraction of element i in the atmospheric gas. Furthermore, 

given the o\ 1 s , the species mass fraction can be determined by FEMP. 


We now proceed to normalize the variables. The distance x,y and the 
distance from the axis are normalized by the stagnation point shock radius of 

curvature R , the velocity u by the freestream velocity u , the density p 

s ® 

by the freestream density p , the pressure p by twice the freestream dynamic 

2 00 

pressure p^u^ , the shock curvature k by l/R g , the total enthalpy H as well 

as the static enthalpy h by the total enthalpy immediately behind the shock 

2 

H g , the stream function <|i by p^u^R , the viscosity by its value immediately 

behind the normal shock p (o) , and finally, the radiative flux Q by 

s 

2R 

s 

pu 3 

r oo oo 

From here on, the equations are all written in nondimensional form. 


Let us divide the shock layer into two 
dominated by the atmospheric gas where f > 0. 
dominated by the ablation gas where \jr < 0. We 
variables § and T| , defined as : 

S(x) = x, T](x,y) = -L } = 

3 ft r 

o s 


sublayers, one is the layer 
The other is the layer 
introduce new independent 


- purdy + pvrkdx (j) 


where 


f3o = 1 for > 0 , 3 q = M for \j; < 0 


‘w s 

2 f ( pv) r dr 2 f ( pv) r dr 

J r WWW J w s s 


M = - 


( 8 ) 


2 

r 

w 
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Hence, M is the total injected, mass flow per unit area as a function of r . 

For a thin shock layer we approximate r w - r s • In Ec 3- ( 8 ), ( P v ) w is the local 
injection rate which remains to he specified. (Note that ( pv) will be negative 

V 

according to the present coordinate system. ) 


By definition, the stream function at shock is \|r = q r and at. wall is 

S CL S 

r 

s 

~J* ( P v ) w r s ^- r s • Hence the range of our new variable T| will be 
o 

from 5 to 


It should be noted that by choosing new variables ^ and 1) the 
continuity equation, Eq. (l), is satisfied automatically, and the y-momentum equation 
can be integrated to give the pressure field as 


P = P c + (PgTl-|)k o U o r 


s s s 


(9) 


Let us now define f s — , then, from the x-momentum equation (Eq. (2)), we obtain 

g 

an equation for f as 


f d / 3f 2 V 
2 sTir si r/ 


\ ± 

9f 2 ± + / 2 

2\ l/ 

t ±\ 

)+ Rep 5 

- p i ff ? ■ M f 

-p) + 

4 - *nP 3 )J 


( 10 ) 


3f 


where f s ^ and the superscript of (±) indicates that different P’s should be 
used dependent on whether T] is greater or less than zero. The 3's are given 
as the following: 

+ , dr s a+ S 1 dU s 

for S: 0 r s / d? > ^2 u d? 

s 


r k r dk 

^-d + ^iV ^ 


P_ k r. 

g + _ _3 + _s_E 

e 4 2 u 


r dr 

o+ _ s s 

-rdT 


( 11 ) 
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for n £ 0 


Pi = 


r M 
s 


dr r JU 

• s s dM 

dg 2 dg 


B- !i^5. 

2 ' "s d ? 


. r k , r dk Mr • 

p: = m (p:-i) + mp: -s. ^ + m m 

3 u 2 1 u dg u dg 

s s = s 


g, k r 
a" _ _3 . S S 
P 4 2 u 


•2 + Mr s dM 

P 5 =M P 5 + ^Tdg 

s 


( 12 ) 


Similarly, the energy and specie conservation equations can be written as 

a r puf shi 2 a r ,, i w 2 af 

an [pr an J 2u s an K^'Pr^ an. 


3” 

+ 1 an I - s^5 - pfij « 


= 0 


(13) 


where Pr is the Prandtl number and Re is the Reynolds number 


Re = 


p u R 

^a> ao s 


and 


H = M 
? “ 


The specie conservation equation becomes 

3C 


a_ 

an 


o ac 
2 V 

pf D - 

H an 

ac 


+ S Reg 

c 5 


2 ^ ar • r i c v? j 


(14) 


where C _ = and D is the effective binary diffusion coefficient. An 

v§ dg 

explicit formula for D will be discussed in Section 4. Sc is the Schmit 


S = 


»s<°> 


c p D (0) 
r oo s 


number defined as 



We note that the nonsimilar terms in Eqs. (10), (13) and (l4) are 
+ + + . , 

P^ff^ , and ’ res P ectivel y- Without these terms, Eqs. (10) to 

( 14) become locally similar equations and § appears only as a parameter. 

The concept of the locally nonsimilar approach is to define those nonsimilar terms 

as new variables and derive equations for them. We thus define 


F h p x ff ( 


Ti 


, H « p,H, , c - p:c 


l C v? 


( 15 ) 


Equations for F, H and are derived in the following manner. For F we 

first differentiate Eq. (10) with respect to ? > then multiply it by P~. An 

equation for F results. In this equation, however, another nonsimilar 
+ 

term p^F^ is introduced. Theoretically one can define this term as another new 

variable and derive an equation for it. A set of an infinite number of equations 

will then result. Truncation must therefore be performed in order to obtain a 

closadset of equations. For the sake of simplicity we will simply neglect the 
+ 

term p F in the F equation. We will call this level of approximation the 
two-equation model. The validity of this two-equation model has been demon- 
strated partially in Ref. 7- In the present work more comparisons will be made 
with some available numerical solutions to examine the validity of the 
two-equation model in treating thin shock layer problems. 


Similarly, the equation for H is obtained by differentiating Eq. ( 13 ) 

± + — 

with respect to § , multiplying it by p^ and neglecting the term P^H^; the 

equation for C is obtained by differentiating Eq. (l4) with respect to ? , 

^ i +— 

multiplying it by P^ and neglecting the term P^C.^ . These equations are given 
as the following 




("“Vi 


B± / 55 


Lf 


31 




+ Reg,. 


+ 

SF ± ± , 2 2s +( ^l^E 

21 i h - hhp - p)- 2 m f + -? 




(16) 
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where 


<"> 5 *§f 


P = ^ 

p ? 


5§ d§ 




35 d§ 


— £ik f M F 9H] cr S_ [ / £M-\ r SH 

BTl[Pr 3T1 J aT|[Pr f dTJ J P 1 an |^Prj an 


+ ““sK’/i In [o“< 1 -f;> f2 lri[( ( "‘< 1 -fe ) ) E f!! | 

**Z h H?* M] * 2 “s f? 

+ 21) 4 + KeS 5 211 Si " B i(pHt)j 4 ’ ° 


a f 2 * c vl a [p 2 DF SCv l+ a* — (d 2 D) f — 

ln[ p sir] + aiu f an" J Pl ^ L Cp D) ? *1 


+2ScRe P 5 n £jj- + ScRe ' v 


= o (18) 


In Eq. (17) the further approximation of neglecting — has been made. 

The explicit formula for the radiative flux divergence will be given in 
Section 3- Needless to say/ it will be a strong function of the absorption 

aQ 

coefficients. To evaluate requires the ^-derivative of the absorption co- 
ol 

efficients. Since a multiband model of the continuum and atomic line radiation 

are being treated, it would have been extremely complicated to obtain the 

ao 

^-derivative of the absorption coefficients. By neglecting ^=- , we greatly simplify 

d| 
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the analysis. An examination of this simplification in terms of the effect of 

9Q 

on the radiative wall flux is considered in Section 5. 

S? 

Equations (10), (13), (1*0, (l6), (17) and (l8) provide six equations 
for f, H, C v F, H, , six unknowns,. These equations can be considered as 
ordinary differential equations with % as a parameter and their solutions can be sought 
locally. They are fully coupled with the radiation field, the mixture thermo- 
dynamics and the mixture transport properties. Methods of computing the radiation 
field and the thermodynamics will be described in the next two sections. 


2.2 Boundary Conditions 

The boundary conditions for f, H, C^, F, H and should be given both 
at wall and at shock. By the definition of f and F , it is clear that at 
shock f = 1 and F = 0. At the wall, a no slip condition holds for the high 
Reynolds number cases of interest. Hence, we find that f = 0 and F = 0. We 
will assume that the shock itself is not radiating and the precursor radiation 
is negligible, then at shock H = 1 and H = 0. At the wall the enthalpy is assumed 
to be given. Thus H = and H = H . Again, due to the high Reynolds numbers 
of interest, the diffusion layer will be much thinner than the shock layer. 

Unless the injection rate is extremely high such that the viscous layer is 
pushed outward near the shock, one would expect that no ablation gas will reach 
the shock. Assuming the injection rate ( pv) is order of p u , then at 

CO 00 

shock C = 0 and C =0. 
v v 


At the wall, the diffusion current for the mass fraction of the atmospheric 

gas must exactly equal the convection current of the mass fraction of the 

atmospheric gas created by the surface mass injection. Hence, at wall we have ■ 

the relation (in dimensionless form) pD & = ( pv) C ReSc where C is the mass 
v K dy K w a a 

fraction of atmospheric gas. Note that .= 1-C^. We have, in the T1 coordinate, 
the boundary condition at wall for as 


0 BC Mr ReSc 
f p D ^ = — (pv) (1-C ) 


(19) 
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The condition at wall for C is found by simply differentiating (19) by 5 then 
multiplying by . We obtain 


o dC ReScM( pv) r 

f < 3 D u ~ + p ? eSc 


^ ( P v) w r s j. 


' e (i-c ) 

e- ' y ' 


|p 2 D + f3^(p 2 D) § 


dc 

arf 


(20) 


2 . 3 Method of Solution 

In attempting to solve Eqs. (l0)-(l8), it was found convenient to trans- 
form the independent variable from T) to a new coordinate 

13 = u, s fptt {21) 

where oo is assigned the value zero at the wall and unity at the shock. 


From Eq. (21) we obtain 


therefore 


T] = <u J fppAu - 

o 


1 

E 


U) 


s 


1 

= -J 

f fpp^u) 

o 


( 22 ) 

(23) 


We define the point uj = (U Q as the value of u) at T] = 0. Then the atmospheric 
gas layer as defined before extends from id = » to a = 1. The injected gas 
layer extends from 00 = ou o to u) = 0. The integration of Eqs. (l0)-(l8) was 
performed in two different directions from oo q toward 1 (3 +, s are to be used) 
and from to to 0 (p 's are to be used). The conditions at the junction 
<D = u) Q were determined by requiring that the values of f, F, H, H, C^, C v 
and their slope in the physical plane be continuous. Hence we require 


df dF dH dOjj. PH dC- 

dy' dy' dy' dy ’ dy ana dy 


be continuous at n = * 

o 


After Eqs. (10) to (l8) were transformed with to as the independent 
Variable and the junction conditions applied, the following formal solutions 
were obtained. 
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For velocity: 


f~~ = G^(<d) + CfG~(u)) + f ((D Q ) 


f = f when (u < u> < 1 , f = f when 0 2 cu ^ uj 
o o 


with 


(D 

G“(cu) = Reg~</ J du> i 


uu 

-2Re3tcu / T|du)" 

, 5 Su,. 




o o 


F+ 4( f2 - 1)- 


-2Repc(u J'lytu)’ 
e «>o do: 


= f 

*0 

[1-G + (1)]G‘(0)M + Gg( l)G”(0) 


f(u>J 


MGg(O) - Gg( l) 


+ l-f(o) Q ) - 0+(l) 


C f = 


G+(l) 


c f = MC f 


For enthalpy profile : N 

H* = g|d)-G^(u>) + CgG^(oj) + H(a> o ) 


H = H when tu^ujSl , H=H when 0 s i m 
o o 


with 


U) ( 1 ) 

o o 



<<») - s %f ( pr - 1 ) f i 4 ” 



;<■> */ 


Gl(to) =| Pr e 

5 


-2ReP u) J PrT)dco 


d<ju 


ub 


H(tu o ) = 


G^(l)[H r - G;(0) + G^(o)] - G’toJCl-G^l) + gJ( 1)]M 


w 3 V 


G*(l) - G~(0)M 


1-G*(l) + Gj(l) - H(cd q ) 
Ch G^(l) 


C H = ” C H 


For specie profiles 

c; = g|(u)) + ctg*U) + c v («» o ) 


(26) 


with 


C = C when m s m s 1 , C=C when 0 ^ u> s <o 

v v o ' v v o 


(0 


r ± 2 
G 6 = ®s 


+ r a c - 2 ScRe e> s X w 

ScRe ^/ W e P ffliC^u,' 


0J n C0_ 


f pd 


„ - 2scite p;j few 

=- e 'u> K 


duo 


< 0 _ 


and 


C ( a) ) 
v o 


9G/-\ dG \ a) ^ MReScr r , l 

>7 ^T/ - *°6 *r] - p5" ~ u ”( pV )w[ G 7 ( 1) _( V l)Gg(0)4MG 6 ( !)G^( o)] 


M 


dG„\ to MReScr r , 

Sr) + W (pv) v [g-(o)m- 07 (i)] 


g 6 (1) + C v* l “o ) 

0*(l) 


C = MC 
c c 
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Similar expressions can be obtained for H, F and as follows 


F = Gg(u)) + G^(o)) + C^G 10 (a>) + F(ud 0 ) 


F = F when tu^^yj^l , F=F when 0 i jj s ^ 


Gg(ou) = - (3* ff 

(On 


( PM-) 


% Sf ,, 

- t— dtu 


PM, Sod 


od 

G^(to) fi: 


+ <P 


° r ± ± 


lS 


-2a, Rep C Tldu)| + + (PM-)- Sf 

5 (o- {2o,„Rep-p-Tl g 

?£%{<* -pY 2 4 ( r < 


i H 55 f 


p 

n 2 „ + 

— + (ju s PixRep 
Sou J 




dOD' 


+ f -2oi ReP s 5 ■‘HdO) 

G 1q =J f e <% cko 


to. 


F(tu o ) = 


G^ 0 (0)M[Gq( 1) + G+(l)] - G^ 0 (1)[G'(0) + G‘(0)] 


G lo^ - G io< o)l51 


+ Gg(l) + G*(l) + F(® 0 ) 


C F = - 


G lo^) 


c f = mc f 


H - = G£q(u)) + G: p2( UJ ) + + H(u> 0 ) 


H = H when ud o £ od ^ 1 , H = H” when 0 s jj ^ ^ 



with 


J 11 


-f I-^u^du) 


t / X F 5H x „ 2 , d -Jf af . SF\ ^ Q f/puA Pr 5H 

X l^ = f 2 au> U s Pr_1 \f d<o au>) PM , do, 


+ 4u u f(Pr-l) + 2u 2 f — /pp,(l-i-jj 

s s? ou) s PM,\^ r Pr 7§ 9(1) 


G“ 2 (cd) = f Prdou / 


+ 013 

-2Re3^u)^y PrT|<iu)" 


0 ) 


I 2 (u>' )do)' 


(IT cu 

o o 


I 2 (U)') = 2Rep^u, s 71I 1 (a)0-2Rep^^u) s 'n 


i tRe6 k f( '“[ pl 5s H+ (^f) l '‘] 


+ r -2ReP ou fTIPrda) 

G 13 = / Pr e P iu o Gu) 


U) 


H(® 0 ) = 


g+ 3 (i)[h w - g^ 1 (o)-g^ 2 (o)] + g‘ 0 (o)[g^ 1 (i)+g^ 2 (i)]m 


g^ 3 (i) - g‘ 3 (o)m 


C H = c ^ 


+ G^(l) + G^(l) + H(d) o ) 

Cg 


Finally 


CT = G~^(<d) + g^ 5 (u>) + c| g i6 (uj) + Cjtuj 


+ + 


V O' 


( 29 ) 


C = C when a s * s 1 . C =C when 0 £ u) ^ u; 

v v o ' v v o 
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with 




SC 

- dm 

00) 


G U */ pD */ 


+ ^ 

-2ReScp-u) s J ^ Tldou" 
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ReScu) 
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e i B 5?l 21 ' 5T - » ^ 


± IF 

2 V(? + P 
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Suj 


dm ' 
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The formal solutions given by Eqs. (24)-(29) are in fact integral 
equations. All the G's are integrals and their integrands are functions of the 
unknowns f, H, C , F, H, C ; the thermodynamics and the radiative flux diver- 
gence Q . The advantages of writing the governing equations into the form 
of Eqs. (24)-(29) are the following: First, the boundary conditions 

are satisfied exactly by Eqs. (24)-(29). Second, Eqs, (24) to (29) can be 
solved by iteration. The iteration scheme is more stable than the integration 
of a second order differential equation with two point boundary conditions. 
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To solve Eqs. (24)-(29) we first guess profiles for f, H, C , F, H 
and . Then, from the thermodynamics subroutine, we obtain the thermodynamic 
properties as well as the transport properties: density, temperature, species 

concentration, viscosity, Prandtl number, and diffusion coefficient. With 
this complete thermodynamic information, the radiative transport subroutine 
will provide the radiative flux divergence Q . • Also, once f, p, p, are known, 

<B , a> , as well as 7] can be determined. Consequently, all the G functions 
can be evaluated and a new (i.e., calculated) value for f, H, C^, F, H, C 
obtained. These calculated profiles will then be compared with the guessed 

ones. If the profiles do hot agree within a specified convergence criterion, 
new profiles will be chosen and the iteration continued until the profiles 
hopefully converge to a solution. 
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3-0 


RADIATION TRANSPORT 


The radiation transport model employed in this study is that of 
Wilson (Refs. 4 and 5). In those references, both the continuum and atomic 
line transport are considered for a system consisting of H, He C, N, 0 atoms 
and significant molecular species. Since helium has high energy level, 
few helium atoms will be excited in the present shock layer problem; their 
contribution to the total radiative flux will be small. The presence of 
helium, however, will affect the composition of gas mixture. Hence, the 
radiation field will be affected by the presence of helium only through the 
the rmodynami c s . 

The detailed analysis of Wilson's radiation transport model can be 
found in Ref. 4. We will not reproduce it here, but will only briefly outline 
the model below. 

It is the radiative flux and the flux divergence that are of interest 
to us. By definition, they are given as: 

00 

i -Si , l(y, v,Q)cos 9dfJdv ( 30 ) 

o Q=4tt 

00 

Q = y.F r =/f li(y,v)[l(y,v,C3) - B(y,v)]dndv (3l) 

o Q=4 tt 

Where l(y, v,fl) is the monochromatic intensity at frequency v and in the 
direction Q , B(y, v) is the Planck function. For the present thin shock layer 
problem, the radiation field is assumed to be one-dimensional. The intensity 
l(y>v,0) can therefore be determined from the one-dimensional radiative transfer 
equations. They are given as: (with the help of the exponential approximation 

+ / -j&/p,(v,y")cLy" 

I (v,y) =J p,(v,y)B(v,y')e y' dy' 

o 
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Where I + is the intensity in the positive y direction and I in the negative 
y direction. A is the dimensional shock layer thickness, £ is the reciprocal 
direction cosine, and p(v,y) the absorption coefficient. With a given thermo- 
dynamic field and p(v,y), the intensity field can be determined. 


The absorption coefficient is separated into continuum and line 
contributions. 

n(v,y) = \i ( v,y) + p, ( v,y) 

The flux may then be expressed as the sum of a contribution due to continuum 
c 

only process F^ and a contribution due to lines but corrected by continuum 
attenuation. 

F = F C + F L (32) 

r r r 


with 


E(0,y) 


E(y,A) 


= r\J" |y*B(v,y' )dE y (y,y' ) -J B( v,y' )dE(y,y’ ) j dv 
00 o 

where the emissive function E_^(y,y') is defined as 

y 

E y (y,y') =1‘- exp[£y* p, C ( v, y")dy"] 

y' 


and 


W,(0,y) 


W,(y, A) 


( 1 

F^ = tt ^ ] f B (y')dW (y'y) -f B (y')dW(yy 

all lines (~ J 


( 33 ) 


( 34 ) 


where the equivalent width W (y y') is defined as 


w ± (yy ' ) = 



vy")dy")l dv 


e 


Similarly, we can write the flux divergence Q into two jiarts 


Q = yF p = Q C + Q L 


( 35 ) 


where 


Q is the energy emitted and absprbed by the continuum 


E v (o, y ) 


r ir" f (y>A) ) 

Q c = 2TT j Ay)\J B(v,y’)dE v (yy') +J B(v,y')dE v (yy') - 2B(v,y)|dv (36) 

o 'o ' 

and 

is the energy emitted and absorbed by the lines. is computed by numeri- 
cal differentiation of the line flux F L . 

- r 


The remaining task is to determine the continuum absorption 
coefficient p C (v,y) and the line absorption coefficient |j, L (v,y) for- 
a particular gas mixture. In Ref. 5, twenty-six-band-grey absorption coefficient 
is modeled for the continuum, and eleven "effective lines" are treated for each 
atom. The detailed frequency range and their expressions for the continuum 
absorption and the definition of "effective line" as well as their frequency 
ranges will not be given here. They are listed in Refs, h and '5. 
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4.0 


THERMODYNAMIC AND TRANSPORT PROPERTIES 


4.1 Thermodynamic Properties 


Thermodynamic variables such as density, temperature and species 
concentrations are needed at each point of the shock layer to solve the conser- 
vation equations and to calculate the radiative flux.. Hence, a relation between 
these variables and the pressure, enthalpy and elemental mass fraction is 
required. For a reacting gas mixture, simple state equations are not available. 

A complete thermochemical equilibrium calculation is necessary to provide the 
needed thermodynamic properties. A free energy equilibrium program (FEMP) 
is therefore incorporated in our computer code ' as a subroutine to provide 
thermodynamic information. The analysis underlying this program is described 
in Ref. 9> The original version of FEMP as described in Ref. 9 does not give 
thermodynamic derivatives. In the present study, certain thermodynamic de- 
rivatives are required to provide the ^-derivatives of the thermodynamic variables 
involved in the governing equations. For example, appeared in Eq. (l6). 
Applying the chain rule, we can write this derivative as 


p § ' a§ dp as 


ah d§ 


— a®. 

t-j'bot. d§ 


(37) 


where c* c is the mass fraction of elemental gas i . 

Equation (37) indicates that we need , |£) and |£ ) 

ap/h,ai dh/p,ai daj/p,h 

We have extended the original FEMP code so that, upon input of p,h and elemental 
gas mass fraction cr , we obtain not only the thermodynamic variables but also 
the thermodynamic derivatives. 


4 . 2 Transport Properties 

Similar to Ref. 4, simplified transport equations of Ref. 10 and 11 
have also been used in this study. The approximate viscosity of a gas mixture 
is given by 


where is 1110 l e fraction and jj,_. is the viscosity of the species i 


and 


^ ~ 2j2 


(■•r 


1 + 


faAi / *jAi 

\ M i/ 


The total thermal conductivity can he written as 


X.k. SC. 

“t E>i ST 

i j J ^-<3 i 


(39) 


where D is the effective binary diffusion coefficient and h is the enthalpy 
of the species i . 


The specific heat at constant pressure of the mixture is 
SC.\ 

C P *E ( V h i §r) ^ 

where C . is the frozen specific heat of each species. For air, Blake (Ref. 13 ) 

pl 

has shown that those approximate equations are valid up to temperatures of 
10, 000° K. 

The temperature dependent diffusion coefficient D used for the 
C, H, N, 0 system was that selected by Hoshizaki and Lasher (Ref. l). This 
curve-fitted diffusion coefficient may not be adequate when helium is included 
in the system. The effect of this uncertainty on the radiative wall flux has 
not been evaluated. 
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5.0 


RESULTS 


5 . 1 Comparisons with Other Solutions 

In order to examine the validity of present locally nonsimilar solutions 
for tnin snock layers, we will make comparisons with other available solutions. 
Unfortunately there is no solution with massive blowing available for comparison. 
The important question of the validity of the locally nonsimilar solution for 
shock layer with mass injection is therefore yet to be answered. 

The following comparisons are made: 

(a) Adiabatic inviscid flow behind a 45° blunted conic shock is compared 
with Olstad's solution (Ref. 6). 

(b) Adiabatic inviscid flow behind a parabolic shock at high Mach 
number is compared with a solution by Schneider (Ref. 14), and 

(c) Viscous, nonradiating flow over a paraboloid at high Mach and high 
Reynolds number is compared with the boundary layer solution by Davis and 
Flugge-Lutz (Ref. 15). All these comparisons are based on perfect gas. The 
results of the comparisons are. present in Fig. 2-5- In Fig. 2, the tangential 
velocity profiles across the shock layer behind the 45° conic shock are shown 

at various locations downstream. In Fig. 3> the body shapes behind the parabolic 
shock are shown. In Fig. 4, the tangential velocity, pressure-, as well as 
the density profiles at z = 10, are shown. Figure 5 shows the convective heat 
flux and the skin friction distribution around the paraboloid. From these 
comparisons it seems that the validity of the present locally nonsimilar solution 
for nonradiating nonblowing thin shock layers has been sufficiently demonstrated. 


5.2 Viscous Shock Layer with Blowing- -but without Radiation 

We present here a typical solution for viscous nonradiating flow behind 
a blunted 45° conic shock with variable surface mass injection. The purpose 
of presenting this solution is to examine the change of velocity and enthalpy 
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Fig. 2. Inviscid Adiabatic Velocity Profiles for Shock Layer 
Behind a 45° Conical Shock 
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Fig. 3. Inviscid Adiabatic Body Shape Behind 
Parabolic Shock at Mach a> 
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Fig. k. Profiles Behind a Parabolic Shock 
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profiles as the flow flows downstream. 


Thermodynamic as well as transport properties are calculated by the 
extended FEMP routine in this solution. We list the necessary input quantities 
as follows : 


Flight Conditions 
Shock Shape 
Shock Nose Radius 
Reynolds Number 
Wall Enthalpy 


= 4.21 x 10 ^ gr/cm^ , u o = 1.6 x 10^ cm/sec 

2 

r 

2 = § 

2 + 1.19 r 

s 

R = 234 cm 
s 

Re = 2.48 x 10 5 

H = 0.2, H = 0 

w w 


Surface Mass Injection Rate 


( p v L = 


w 


1 + r 


Atmospherical Gas : Air Injected Gas : Air 


(4l) 


The results are plotted in Figs. 6 and 7- Figure 6 shows the shock and body 
shapes. Also plotted in the figure is the interface between the atmospheric gas and 
the injected gas. This interface is defined as the surface where the stream 
function is zero. We see that the body shape is not significantly different from the 
shock shape. The injected gas layer grows but at a slower pace than the growth 
of the total shock layer. This is because of the decreasing surface, injection 
rate around the body. Figure 7 shows the tangential velocity and static enthalpy 
profiles. It is seen that the viscous layer is detached from the wall due to 
mass injection. Hence, both the convective heat transfer and the skin friction 
are practically zero. The inviscid layer near the wall becomes thinner as the 
flow flows far downstream. 


5-3 


Inviscid Shock Layer with Radiation but without Blowing — Evaluating the 
Effect _ 

In the derivation of the equation for H we have neglected the term 


k 
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Fig. 6. Shock and Body Shapes for a Blowing Case 
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Fig. 7. Enthalpy and Velocity Profiles for a Blowing Case 
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90 

, namely, we have neglected the contribution of the streamwise rate of change 
9s . 

of radiative flux divergence to the energy equation. In order to assess the 

90 

effect of ^ on the radiative wall flux we have made computations with and 


without this term. The computation is for an inviscid, nonblowing case; the 

flight condition and shock shape are the same as given by Eq. (4l) and the 

shock nose radius is 100 cm. The three-band radiation transport model for air 

continuum as described in Ref. 7 is employed for this particular case. Figure 

8 shows the result for radiative wall flux, the dotted line is the result for 
3 Q 30/ 

^ = 0, while the solid line represents the result for ^ f 0. We see that the 
90 do . 

result for = 0 is always higher than for f 0 . The maximum difference 
o§ 9§ 

between these two results in this case is about 15$. A 15$ inaccuracy seems to 

be reasonable to accept in view of the simplification it will introduce by 
90 £0 

neglecting this term. To evaluate ^ in the present study would require the 
d S 9| 

^-derivatives of Eq. (35)- It involves, first, the transformation of those 

equations from (x,y) plane to (§,T|) plane, then taking partial derivatives with 

| by keeping "q fixed. A tremendous amount of algebraic manipulation would 

have to be made in this process, and this can all be avoided by simply neglecting 
90 — dQ 

— in the H equation. Although the effect of on the wall flux is not 
negligible, we have demonstrated for the case of three-band radiation models 
that the effect is on the conservative side and is reasonably small. 


5 . k Completely-Coupled Radiating Shock Layer with Blowing 

A sample calculation is made for the shock layer problem of a Jupiter 
probe. The flight conditions, as well as other necessary input data, are 
listed below: 

Flight Conditions x I® gr/cm , = b x 10 cm/sec 

2 

r 

g 

Shock Shape Z = g + 1.88 r 

s 


Shock Nose Radius 


R = 20 cm 
s 


Wall Temperature 


T = 4500 °K 
w 
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-0.4 


Surface Mass Injection Rate (pv) = 

w 

" ' ' 's 

Atmospherical Gas: 7^ Hydrogen, 26$ Helium; by mass. Injected Gas: Carbon 

Radiative Model: Given by Refs. 4 and 5- 


Based on the flight conditions listed above, the Reynolds number as 
ilculated by FEMP is 9-267 x 10 ^ and the density ratio ( p o0 /p s ) across the 
Drmal shock is 0.0964. 


In this calculation, the ^-derivatives of the thermodynamic transport 

uch as and ^ are neglected for the following two reasons : 

oz oz oz 

(a) Those derivatives turn out to be very sensitive to the local gas 
omposition [hue to the very steep gradient of the ablating product concentration 
rofile near the interface (the viscous layer)]. This sensitivity causes the 
teration to become unstable. 

(b) The effect of and should be confined within the viscous 

oz oz 

ayer; in the present high Reynolds number cases this layer is indeed very 

•hin. Hence, the effect of neglecting ^ , PP and ^ on the wall radiative 

0 § 0 § 0 § 

.'lux should be small. 


In order to examine the effect of ablation on the radiative wall flux, we also 
lake a calculation for the nonblowing case. The results are shown in Figs. 9 - 
14. 


In Fig. 9, the shock shape, body shape, as well as the interface, are 
shown. The interface is defined as the surface where the stream function 
Ls zero. It is found that I^/R^ — 0.42. In the present case, R s = 20 cm. 
Hence, the body nose radius of curvature R^ is approximately 8.2 cm. The 
asymptotic body half cone angle is approximately 57° • Also shown in Fig. 9 
is the body shape when there is no blowing. This body shape is practically 
identical to the interface shape for the blowing case. The last point in the 
present calculation is at r^/R^ — 2, r ^/^ s — 1-85- At this point we see that 
the shock angle is close to its asymptotic angle which is 62° . Also, we see 
that, at this point, the shock layer is not exactly a thin layer. 
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Fig. 9- 


Shock and Body Shapes for a Radiative Coupled Viscous 
Shock Layer with Blowing 




Figure 10 shows the surface radiative flux for both blowing and no 
blowing cases. -A drastic reduction in radiative flux due to mass injection is 
seen. The flux distribution for the no blowing case shows an anomalous behavior, 
namely, the maximum heating is not at the stagnation point. For the present 
particular case, the downstream flux maintains at a higher level than that 
at the stagnation point. In order to explain this behavior, we have plotted 
the shock layer thickness distribution and the temperature profiles in Fig. 11. 
The shock layer thickness increases rapidly near the nose followed by an almost 
constant rate downstream. By comparing the temperature profiles at Z = 0 
and Z = 0.03, we see that the temperature changes very little. On the other 
hand, the shock layer thickness increases by about 40 $ from 1.2 to 1.8 cm. This 
results in a higher flux at Z = 0.03 than that at Z = 0. The temperature 
decreases as Z increases, but remains at a high level as one can see from 

the profiles. The decreases in shock layer temperature will result in a lower 
surface flux, but it is compensated by a thicker shock layer. Hence, the net 
result is the flux maintains at a nearly constant level which is higher than 
the stagnation value. For blowing cases, the flux decreases monotonically 
as can be seen from Fig. 10. The heating level for a blowing case is mostly 
controlled by the injected gas layer. A monotonically increasing injected 
gas layer should result in a monotonically decreasing flux distribution. 

Figure 12 shows the flux reduction due to the mass injection. The reduction 
ranges from about 40$ at stagnation region to about 70$ far downstream. This 
indicates that injected carbon gas is an effective heat shield in the sense 
that most of the shock layer hot gas is dumped into the wake. A detailed 
examination of the spectral distribution of surface radiative flux reveals that 
most of the continuum flux reduction occurs in the carbon molecular bands of 


and C^. The reduction in line flux is also a significant portion of the 
total flux reduction. It is interesting to note that from Fig. 10 one sees 
that at the stagnation point, the line reduction accounts for about 60 $ of the 
total flux reduction. 


Figure 13 shows the typical profiles for tangential velocity, enthalpy 

and injected gas mass fraction profiles at r = 1.22. It is seen that about 35$ 

s 

of the shock layer is an injected gas layer. The shear layer is thin, about 10$ of 
the shock layer. Most of the injected gas layer is inviscid. The temperature 
profile and the C, C^, concentration profiles are shown in Fig. l4. 
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COMPUTED POINTS 
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Fig. 11. Shock Layer Thickness and Temperature Profiles for 
the No Blowing Case 
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Fig. 12. Surface Flux Reduction Due to Blowing 
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Fig. 13. Temperature, Velocity and Injected Gas Mass 
Fraction Profiles at Z = 0.35 
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Enthalpy, C^, and Profiles at Z = 0-35 
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Fig. 14. 


6.0 


CONCLUDING REMARKS 


Viscous radiating shock layer flow over a blunt body with massive 
blowing has been solved by a locally nonsimilar approach. Detailed equilibrium 
computations for the thermodynamic and transport properties are included in 
the solution. Molecular band and atomic line, as well as continuum radiation, 
are considered. 

Comparisons with prior inviscid solutions for flow behind blunted 
conic shocks and flow behind parabolic shocks show the present solution method 
adequately treats nonsimilar aspects of flow around blunt bodies. A sample 
calculation for a typical Jovian entry problem has been made. A ho $ reduction 
in total radiative flux is found at stagnation point due to the injection of 
carbon gas. The reduction increases to about 70$ at about four body nose 
radius downstream. This increase of flux reduction indicates that most of the 
radiative energy is absorbed by the injected carbon gas and dumped into the 
wake . 


An examination of the spectral distribution of the radiative flux 
reveals that most of the flux reduction occurs in molecular bands of C ^ and 
C^. Due to the large Reynolds number and massive local blowing rate, the 
solutions indicate that the shear layer remains thin (about 10$ of the shock 
layer). The injected gas layer is mostly inviscid. 


Due to the very steep gradient of the ablating product concentration 
profile near the interface (the viscous layer), the ^-derivatives of the thermo- 
dynamic transport such as 


^ an( i are very sensitive to the local gas 


. as ' a? as 

composition. This sensitivity causes the iteration in the numerical computation 
to become unstable. Therefore, in the sample calculation for the Jovian entry 
problem, S-dnrivatives of the thermodynamic -transport have been neglected. In 
the present higher Reynolds number and massive blowing case, the effect of 
neglecting those derivatives on the wall radiative flux should be small. 
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